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Abstract 

This paper is a contribution to the development of a framework, to be used in the context of 
semiclassical canonical quantum gravity, in which to frame questions about the correspondence 
between discrete spacetime structures at "quantum scales" and continuum, classical geometries at 
large scales. Such a correspondence can be meaningfully established when one has a "semiclassical" 
state in the underlying quantum gravity theory, and the uncertainties in the correspondence arise 
both from quantum fluctuations in this state and from the kinematical procedure of matching a 
smooth geometry to a discrete one. We focus on the latter type of uncertainty, and suggest the use of 
statistical geometry as a way to quantify it. With a cell complex as an example of discrete structure, 
we discuss how to construct quantities that define a smooth geometry, and how to estimate the 
associated uncertainties. We also comment briefly on how to combine our results with uncertainties 
in the underlying quantum state, and on their use when considering phenomenological aspects of 
quantum gravity. 
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I. INTRODUCTION 



Ever since quantum mechanics and general relativity became two of the pillars of contempo- 
rary physics, it has been clear that a new, more general theory is needed that contains both 
of them in some appropriate limit. The task of finding such a quantum theory of gravity has 
indeed been an excruciating one and has not yet been completed, in part because of the nu- 
merous conceptual challenges that such a theory faces [ij. One of them is the question that 
will motivate our considerations: What is the fundamental structure of spacetime? That 
is, can we expect the continuum picture of a (smooth) manifold to be valid at arbitrarily 
small scales? Do we need to modify it into a "foamy" picture with a manifold of fluctuating 
topology and possibly dimensionality? Or should it be entirely replaced by a different, 
quantum picture that is discrete, polymeric, fuzzy? Is it even appropriate to pose the basic 
questions in such a geometric language? It is by now generally believed that the smooth 
manifold picture is inadequate, but there is a wide variety of approaches to the way in which 
the quantum realm manifests itself, which in some cases depend on the prejudices of individ- 
ual quantum gravity practitioners, for example regarding the amount of classical structures 
and non-dynamical elements that remain in the final formulation. In the present work we 
shall not pretend to be free from such prejudices; however, we will attempt to attack our 
main problem from a general perspective, and present a framework that could be adapted 
to different approaches and formulations of the dynamics. 

Within the search for a consistent theory of quantum gravity, there has been over the past 
several years an increase in the amount of work devoted to the semiclassical sector of the 
theory and its physical predictions. These subjects involve several issues, such as identifying 
those states which have a semiclassical interpretation, and predicting phenomenological ef- 
fects due to these states that can be seen as quantum gravity corrections to the dynamics of 
quantum fields. As a result of their different uses of the classical structures, each approach 
faces a different challenge when it comes to the question of the (semi-)classical limit. As loop 
quantum gravity is the main inspiration for our framework, let us look at the semiclassical 
limit there. 

A key aspect, underlying the very possibility of phenomenological predictions, is that of 
establishing a correspondence between the structure in terms of which quantum gravity is 
formulated, and the smooth geometry in terms of which the usual field theories are formulated 
and observations interpreted. One expects the two descriptions to be very close at large 
scales, and to start departing significantly at small length scales, since all of the field theories 
we use to describe the behavior of matter, and even of gravitation, are based on manifolds 
that become "featureless" at small scales, and are actually locally fiat in an infinitesimal 
neighborhood of each point. We will therefore take the point of view that the manifold Ai 
is only a convenient tool which gives us an effective "low-energy" picture, within a certain 
degree of accuracy at each scale. With this observation at the center of our discussion, we 
will describe a framework that will quantify this accuracy, and identify scales at which one 
can cross over between the quantum picture of spacetime, assumed to be discrete, and a 
continuum, classical geometry. 

The issue of recovering a continuum theory as an appropriate limit of a discrete one is of 
course not new, and has been studied in various contexts for a long time. In quantum gravity 
the best-known approaches in which the issue arises are the ones based on piecewise linear 
manifolds, such as the various versions of Regge calculus Q or dynamical triangulations 
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graph-based approaches such as loop quantum gravity j3[, the more recently developed spin 
foam models or causal sets 0, but the list includes many others, both covariant and 
based on space+time splittings. In the causal set approach and in dynamical triangulations, 
the only variable used is combinatorial, while in other approaches one uses a combinatorial 
structure which is "dressed" with additional variables, such as edge lengths of a manifold 
triangulation in Regge calculus, or holonomies of connections along graph edges in loop 
quantum gravity. Although there may or may not be a "minimal length" in those frameworks, 
one may think of their discreteness as giving rise to a characteristic length scale I (whose 
value in terms of fundamental constants is to be determined within the theory) above which 
the discrete structure VL can be considered as kinematically well approximated by a smooth 
manifold. 

When one faces the issue of recovering the continuum theory, the first question that 
arises is whether the discrete dynamics converges to that of the continuum as £ 0. To 
address this question, one usually considers a fixed smooth manifold with metric {M.,gab) 
(to be thought of as either the spatial or the spacetime geometry), and embeds in it an 
unspecified, but increasingly finer sequence of discrete structures VL{t) (for example, trian- 
gulations of M. with edge lengths induced by Qah of order £), and shows that the value of 
the discrete version of some quantity Q that governs the dynamics, such as the action in a 
covariant approach, approaches the value of the continuum version, calculated for {M.,gab), 
as Q(f2(£)) = Q{A4,gab) + C(£^). A number of results of this type are known [Bj. For the 
most part, the emphasis has been on features of the discrete dynamics itself and on showing 
that the discrete variables approach the continuum ones fast enough. 

Here we want to explore the kinematical aspects of this correspondence in more detail. As 
already mentioned, we feel that, given the importance of potentially observable corrections to 
continuum theories, it is important to have a way to quantify the extent to which a discrete 
geometry corresponds to a smooth one, and the amount of uncertainty in the correspondence. 
In other words, if the discrete theory is fundamental and not an approximation, the limit 
i —>■ is not to be taken, and we need to know exactly what the 0{P) terms are. This is the 
goal of the framework that we propose. In this work, the first in a series of papers, we intend 
to motivate the issues to be addressed and to provide the first steps. We start by considering 
only those aspects that are related to an underlying combinatorial structure; aspects related 
to additional variables to be assigned to the underlying discrete structure will be treated 
separately in a forthcoming publication |^. All considerations will be independent of the 
dynamics of the theory, and in fact we don't need any details of the quantum theory, except 
for occasional references to an underlying quantum state assumed to be semiclassical. 
This paper can be separated in two parts. In the first one we develop a way of assigning a 
smooth classical geometry to a given (random) graph. The second part places this procedure 
in the more general context of building a macroscopic geometry {Ai,g), starting from a 
semiclassical state "^sc, or a graph Q together with the expectation values of observables 
O, and raises the related issues of coarse-graining and combining statistical and quantum 
uncertainties. 

Let us now comment on the way in which our work differs from previous results on 
the continuum limit of a discrete theory. The first key feature of our approach is related 
to the fact that, rather than reconstructing a given smooth geometry, we are interested in 
constructing one using a discrete structure Q, which may not even be embeddable in a smooth 
geometry. For example, in the picture of quantum geometry arising from loop quantum 
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gravity, the geometry at Planck length appears to be distributional, with support on the 
edges of a graph. In the current formulation of this approach, those graphs are embedded 
in a given manifold A4, but the expectation is that it should be possible to formulate it 
in terms of abstract graphs. The theory should then include criteria for recognizing graphs 
which look like manifolds at large scales, and specify how to determine the emergent geometry 
and evaluate the uncertainties involved in the construction. Several scales, associated with 
qualitatively different descriptions of the geometry, will arise in this process, and it should 
ultimately be possible to establish a correspondence between those provided by continuum 
theories, such as the Planck length Lp at which classical geometry is expected to break 
down, and the ones provided by the discrete theory, which in the present paper are purely 
combinatorial and dimensionless, such as the amount of coarse-graining necessary for Q to 
be embeddable. 

There is a caveat, however: a random complex carries information from the metric of 
the underlying manifold. As will be shown in detail, curvature quantities can be related 
to combinatorial properties of the complex and vice versa. But this raises a potentially 
worrisome point: would it be consistent to use a complex that determines, in the above 
sense, a classical metric gi, say, but dress it with semi classical states that are peaked around 
a completely different, non gauge-equivalent classical metric 5(2? To circumvent this potential 
problem, it would be convenient to have criteria to decide on the right random complex to 
use for the semiclassical situation at hand. This will imply, at this level, to use the quantum 
geometry g2 implied by a quantum state \l/ and compare it with the geometry gi consistent 
with the discrete structure fl. 

As tools for constructing a geometry {M.,gab), we will identify examples of quantities Qx 
associated with appropriate regions or submanifolds X of Ai; the second key feature of our 
approach concerns the way in which uncertainties AQx are calculated using statistical tech- 
niques. Given one discrete structure Q, the AQx represent the uncertainty in the estimate 
of the effective geometry {A4,gab) that Q could be considered a discretization of. Such sta- 
tistical fluctuations can be, and have been, calculated in some cases, using both analytical 
and numerical methods. Here we will use as far as possible analytical techniques, in the 
spirit of the random lattice approach to gauge theories in Minkowski space pioneered by T D 
Lee and collaborators |^. Most of the explicit calculations will be done in two dimensions; 
results in three or more dimensions in general will have to be obtained numerically. 

More specifically, in SeclHlwe define the setting of our work, introduce and motivate the 
choice we make for Q, that of a cell complex, and comment on the possibility of obtaining 
cell complexes starting from graphs. Sec IIIII contains the main results in this paper; after 
an introduction, clarifying what we mean by constructing an approximate geometry from Q, 
we show explicitly how to carry out the construction in 2D. In Sec II VI we return to more 
basic questions about our discrete structure: we discuss the role of different length scales 
in the discrete-continuum transition, and the more general setting of structures that are 
not embeddable at a certain length scale. For example, discrete structures associated with 
candidate semiclassical states for quantum gravity will probably not be generic graphs, and 
it seems reasonable to start exploring the obstructions to graph embeddability by looking 
at cases that are not too "severe", where the discrete structure is, in an appropriate sense, 
"almost embeddable" . To that end, we discuss the notion of coarse-graining for a cell com- 
plex, intended to produce in those cases an embeddable, smoothed-out version of the discrete 
structure. 
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In Sec El we discuss the consequences of the quantum fluctuations (AQ/)*, and possi- 
ble uses of our results. As stated earlier, from the point of view of the continuum, those 
fluctuations contribute to the total uncertainty in the geometry {J^,gab)', thus, on the one 
hand they will ultimately allow us to quantify the goodness of \E' as a semiclassical state, 
and on the other hand they will play an important role in the relationship between \E' and 
continuum-based phenomenology. Note that these modifications are of a different nature 
than the corrections to field dispersion relations that arise purely from discretization ef- 
fects in Ref Finally, we return to the motivation for this work, and discuss possible 
applications in loop quantum gravity phenomenology and other directions for future work. 

Regarding notation, we will follow the following convention. Statistical averages and 
quantum expectation values will be indicated by angle brackets, as in (Q), and means with 
respect to probability distributions by overbars, as in Q; as for uncertainties, (AQ)|, will 
denote a quantum fluctuation, while {AQY (the subscript "c" being understood) or ctq will 
denote the statistical uncertainty or variance of a classical probability distribution. 

II. MATHEMATICAL SETTING 

Our framework can be seen as a "bridge" between a discrete, pre-geometrical description of 
spacetime, motivated by quantum gravity ideas, and the classical, continuum-based geometri- 
cal description, for situations in which the quantum theory provides us with a "semiclassical 
state". In this section, we define the notions needed to translate this statement and the 
conceptual points discussed in the introduction into a specific program. 

A. Cell Complexes and Tilings 

The most basic variable in this paper will be a cell complex; given their central role in what 
follows, we start by recalling a few useful definitions and facts about cell complexes and 
their relationship with manifolds. In topology, a /c- dimensional (open) cell is a space homeo- 
morphic to the interior of a /c-ball. A cell complex is a set of nonempty, pairwise disjoint 
cells, such that (a) The closure of each cell is homeomorphic to a ball and its boundary to a 
sphere in some dimension, and (b) The boundary of each cell is a union of cells; in our case, 
this will always be a finite union. (0-dimensional balls are pairs of vertices.) 

Given a differentiable manifold A4, a cell decomposition of is a cell complex homeo- 
morphically embedded in it. Our assumptions then imply that the cell decomposition is 
locally finite, in the sense that every compact subset of A4 intersects only a finite number 
of cells. For example, a finite 3-dimensional cell complex Q (one whose maximal cell di- 
mensionality is 3) consists of a set of A^^q vertices vj, Ni edges ej, N2 2-cells uj, and 
3-cells C/; if is a cell decomposition of Ai, then the 3-cells together with their boundaries, 
(U/ Ci) U (U/ uji) U (U/ e/) U (U/ w/) ~ A^. We will say that a cell complex is embeddable 
if there is a differentiable manifold of which is a cell decomposition. 

We will also mention the more general concept of tiling of Ai. This term often denotes 
a collection of pairwise disjoint open subsets of A4 whose closures cover A^; here we 
will consider a tiling to include the union of the boundaries of the uj, partitioned into 
submanifolds of various dimensionalities. A simple example will illustrate the concept. Given 
a smooth loop a in S^, such that S^\q; consists of two open "half-spheres" uji and ^2, the set 
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{tui, to'2, a} is a tiling of S^, but not a cell decomposition. However, if we pick a point v ^ a 
and call e the edge a\v, then Ci := {ui, uj2, e, v} is a cell decomposition of S^; alternatively, 
if we pick two points on a and join them with an extra edge that does not meet a elsewhere, 
we divide into three open wedges which, together with the elements on their boundaries, 
make up another cell decomposition €2- In a cell complex Q homeomorphic to a D-manifold 
Ai (but not in any tiling), the Nic{Q) satisfy 

= ^(^)' (2.1) 

where x(^) is the Euler number of the complex Q, or of the manifold A4. 

However, cell complexes need not be embeddable. They could, for example, have "regions" 
with different dimensionalities; in an embeddable cell complex every cell is, or is on the 
boundary of, one of maximal dimensionality. For a different type of example, consider the 
above cell complexes Ci and C2, and remove their 2-dimensional cells; in the first case the 
remaining 1-dimensional skeleton is homeomorphic to S^, while in the second case one is left 
with a non-embeddable complex. 

A useful operation on cell decompositions is duality, which produces a new D-dimensional 
cell complex Q* from any given Q of the same dimensionality. One associates with each k- 
dimensional cell uj in Q (a D-dimensional cell complex must have cells of all dimensionalities 
< k < D) a — fc)-dimensional dual cell u* in Q*, whose boundary consists of the duals of 
all cells which have u on their boundary. If is a cell decomposition of a manifold Ai, then 
fl* is also homeomorphic to A4; however, since the duality f2 fi* is an operation between 
abstract cell complexes, in general there is no natural embedding of Q* in Ai (the mapping 
f : Q ^ Ai does not induce a mapping /* : Q* ^ Ai, unless At is endowed with more 
structure). Duality is defined for general tilings, but their duals may not be homeomorphic 
to the original Ai, while non-embeddable cell complexes may not have well-defined duals. 



B. Triangulations and the Voronoi Procedure 

Of all types of cell decompositions of manifolds, the most useful ones for us are triangulations, 
in which the cell complexes are simplicial, and their dual complexes. In triangulations, of 
course, all 2-cells (triangles) have 3 edges and 3 vertices, all 3-cells (tetrahedra) have 4 faces, 
6 edges, and 4 vertices; in general, all fc-simplices have k + 1 faces on their boundary, etc. 
Their dual cell decompositions therefore satisfy incidence properties which state that each 
{D — k)-cell is on the boundary of (is shared by) k + 1 cells of dimensionality one unit higher, 
etc, and can be concisely written as follows: For each /-dimensional cell u & Q, the number 
of /c-cells that have u in their closure (with < I < k < D) is 



Nk\i{u^) =( u 1 ]■ (2-2) 




In particular, each vertex has A'^iio = D + 1 edges. Thus, in two dimensions all dual vertices 
are trivalent, while in three dimensions they are shared by four edges. This property already 
makes such complexes useful, since for example quantum geometry results in loop quantum 
gravity show that 4-valent vertices of graphs in 3 dimensions are the fundamental units of 



volume Also, if each edge terminates at two vertices, we find a useful relation between 
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the total numbers of edges and vertices (if finite), 

N,{n) = l{D + l)No{n) . (2.3) 

However, the main reason why these two complexes are useful is more general; it lies in the 
fact that they can be obtained in a manifold using just a set of points as input, and in the 
way they encode geometrical information on (A^, g^b) when the set of points they are based 
on is chosen at random. 

Given any locally finite set of points pi in a Riemannian manifold {A4, Qab), one can obtain 
from it a tiling of Ai; if the points are at generic locations and sufficiently dense (with respect 
to both local length scales, determined by the metric, or global ones, determined by the metric 
and topology), the result is actually a cell decomposition, called the Voronoi complex, and 
its dual is called the Delaunay triangulation. We start by introducing the procedure in 
general, without any additional assumptions on the set of points, and then consider in the 
next subsection the case in which the {pi} are randomly sprinkled points. (Capital latin 
indices /, J, will be used to denote points or elements of various dimensionalities in a 
complex; the type of object they refer to should be clear from the context.) 

For each embedded point pj, we can define an open region uj C as the set of all 
manifold points which are closer to pi, with respect to g^b than to any other pj] clearly 
the union of the closures of such regions is A^, so the {wj} define a tiling. The boundaries 
of the uj are made of manifold points that are equidistant from more than one of the 
p/'s; those equidistant from pi and pj and closer to them than to any other px are the 
codimension-1 common boundary of the D-dimensional regions ui and Uj around pi and pj, 
respectively. Common portions of boundaries among more than two cells, if present, define 
lower- dimensional portions of the dui^s. We will always call the resulting complex a Voronoi 
complex, even when not all the sets just described are cells; if the pi are close enough to 
each other compared to all length scales in the manifold, associated with the metric or the 
topology, we actually obtain a cell complex Q homeomorphic to A4. 

Generically, vertices of a Voronoi complex are equidistant from D + 1 points, since a 
higher number of points at the same distance can only be obtained in degenerate situations; 
then the -D + 1 ways of picking D of those points define D + 1 Voronoi edges incident on 
that point. In fact, when the pi are at generic locations, all of the relationships ()2.2|) are 
satisfied, regardless of whether we have a cell complex; because in our applications the points 
will be randomly sprinkled in A4, we will not discuss special arrangements and will always 
assume that those relations hold. Given any Voronoi vertex vj, the D + 1 points it is close 
to are on the surface of a Z)-sphere around vj and define a simplex u}, dual to vj. If Q 
is a cell complex, there are enough Voronoi vertices for the dual complex Q* defined by all 
those simplices to be a Delaunay triangulation of the manifold, which has the pj themselves 
as vertices. Intuitively, k + 1 sprinkled points that are "clustered closely enough" define a 
fc-simplex in Q* that lies in the /c-plane through those points. Thus, the concept of Delaunay 
triangulation is less general than that of "Voronoi tiling" . 

When does the Voronoi procedure not produce a cell complex? A simple example is that 
of a 2-sphere with two points pi and p2 chosen on it, in which case the procedure gives the 
tiling we called {cui, 0^2, a} in Sec 111 Al and its dual consists of pi, p2 and an edge between 
them, which is not homeomorphic to S^. In this example, both the non-trivial topology 
and the high curvature (in terms of the point density) contribute to the outcome. However, 
one can easily modify it into other examples in which only one of the factors is present 
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FIG. 1: Example of 2D Voronoi (thin lines) and Delaunay (thick lines) complexes; the points they 
are based on are the Delaunay vertices. Notice that in 2D edges are dual to edges, whereas vertices 
and 2-cells are duals of each other, and that dual edges are orthogonal to each other, although they 
do not necessarily meet. 

(a fiat cylinder, or a plane with a small region blown up into a long tube or a balloon, 
respectively), and it produces a similar effect. All such Voronoi tilings satisfy the incidence 
relations (j2.2j) . and those same relations, imposed on an abstract complex, guarantee that 
it is homeomorphic to some manifold. For the time being, however, we would like to work 
with ordinary cell complexes. This means that for us an embeddable complex will be one 
such that the relations ()2.2|) are satisfied and in which every cell has on its boundary (the 
appropriate number of) cells of all lower dimensionalities. 

By construction, the two types of complexes just defined are embeddable. However, 
since in this paper our intention is to construct geometries from discrete structures, we are 
interested in characterizing those structures which can arise from the Voronoi procedure 
in some Riemannian manifold. This condition excludes the "lower-dimensional parts" in 
the complex, but other possible "defects" are not excluded. Voronoi complexes provide the 
most convenient combinatorial set of conditions for being embeddable in a manifold. These 
conditions exclude the "degenerate" Voronoi complexes mentioned earlier (which is not a big 
loss), but they do include ones which are not locally finite, and could easily be extended to 
the generalized complexes that are not made of topological cells. 
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C. Random Voronoi Complexes 



In order to understand how discrete structures embedded in a manifold Ai encode its geom- 
etry, we need to introduce a few notions related to randomly distributed points. A random 
point distribution on Ai with a volume element (in particular, one given by a Riemannian 
metric, ^/g) is the outcome of a uniform, binomial or Poisson, point process. If the total 
volume V_\4 is finite, a uniform point process is specified by stating that, each time a point 
X is chosen in Ai, the probability that x fall in any given measurable region X <0 Ai is 

P(x G X) = Vx/Vm , (2.4) 

or the infinitesimal version, that the probability density is pMi^ly/g) = Vm y/g (^^x. If the 
process is repeated N times, with no correlations among points (we will not keep track of 
the order they came in), we get a uniform sprinkling of points with density p := N/Vm that 
for our purposes we will think of as being of order C.~^ . The probability density for each 
point to fall in an infinitesimal region d^x is 

d/i = py^d^x. (2.5) 

One of the most useful finite probabilities in this context is the one for exactly k points out 
of to fall inside X (without specifying which ones). It is easy to see that this probability 
follows a binomial distribution, 

which, as Vm and become very large, with p = constant, approaches a Poisson distribution, 

Q-pVx (pVx)'' 

P{k,X \ N,M) ^- 



This last equation justifies the name Poisson distribution that is often used for the sets of 
points used in this paper, and corresponds to the infinite volume situation. 

A random Voronoi complex in a manifold {Ai, Qab) is the result (assumed to be a cell 
complex for the time being) of applying the Voronoi procedure of Sec III CI to a Poisson 
distribution of points in {Ai, gab) (see, e.g., Ref and references therein). This is the 
discrete structure we use in this paper to encode information on a spatial geometry. Such 
structures have been called "random lattices" in the context of gauge theory (the subject 
is actually older than that, but for references with physical motivations somewhat related 
to ours, see Refs and ^3). The randomness of the point distributions and their finite 
density imply that all complexes are locally finite, and all vertices are {D + l)-valent with 
probability 1. Most of the general discussion will be valid for complexes and manifolds of 
arbitrary dimension D, but actual calculations will be carried out for 2-dimensional ones. 



D. Remarks on the Use of Voronoi Complexes 

The discrete structures one uses in loop quantum gravity to construct spin networks, the 
basic states that form the usual basis for the kinematical Hilbert space Hkm, are graphs 
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embedded in a manifold. Thus, although we don't know yet how to formulate a manifold- 
independent theory, it may be useful to try to establish a correspondence between certain 
types of graphs and the Voronoi cell complexes that our statistical machinery is based on. 
Thus, suppose that we are given a graph 7, consisting only of vertices and edges, but without 
higher-order cells. We then pose the following question: Can we construct a full Voronoi 
complex just from 7? To recover Q as an abstract complex in D dimensions, what we need 
to do is specify which edges in 7 form (the boundary of) an elementary 2-cell, which of the 
latter form (the boundary of) a 3-cell, and so on. If all goes well, the resulting cell complex 
will satisfy the correct incidence relations for a Voronoi complex, as specified by ()2.2|) . 

Let us begin with a proposal for a construction in two dimensions. Given a graph 7 in 
which every vertex is trivalent, define a loop to be a chain of consecutive edges 6162 ■■■ = 
(f 1 f 2 ^ "f^x ^ f^K+i) that closes on itself, i.e., Vi = vk+i- Most loops are not 

to be thought of as boundaries of 2-cells; we call plaquette a loop a such that, for any two 
vertices vj and f j G a, the shortest path in the graph between vj and vj is part of a. The 
cell complex Q we are looking for has as 0-cells and 1-cells the same ones as 7, trivially, and 
its 2-cell are identified with a set of plaquettes a such that is a Voronoi cell decomposition 
of a 2-manifold, i.e., such that every edge is shared by exactly two 2-cells; if such a choice is 
not possible, we consider the graph non-embeddable. However, we expect that if we apply 
this construction to the 1-skeleton of an actual Voronoi complex Q, we recover the original Q. 
The main questions will then be: How do we recognize from the graph whether it corresponds 
to a "good situation"? What can we say about cases where it does not? 

Similarly, to construct a 3D Voronoi complex from a four-valent graph, we start by defining 
a candidate 3-cell C as a finite set of plaquettes {ai, a2, ■ ■ ■ , such that (i) every edge is 
shared by exactly two 2-cells and aj, as in the 2D construction above, (ii) the topology 
of C is that of a 2-sphere, and (iii) for every two vertices Vi G ai and vj G aj, the shortest 
path in the graph between them is part of C. A collection of 3-cells defined in this way gives 
a good Voronoi complex if all edges in it are shared by exactly three plaquettes, and each 
plaquette by exactly two 3-cells. These definitions could then be generalized in an obvious 
way to higher dimensions, but in this paper we will not need to consider explicitly cells of 
higher dimensionality. We can call abstract D- dimensional Voronoi graph one such that all 
vertices have the same valence D + 1, and the above construction gives a good D-dimensional 
Voronoi complex. 

Let us conclude with two remarks on the concepts we have introduced so far. First, we 
are not suggesting that all semiclassical quantum gravity states are associated with Voronoi 
complexes, just as in ordinary quantum mechanics, not all semiclassical states are coherent 
states. However, the latter have properties that make them easy to work with, and they 
encode in a convenient, minimal set of parameters a point in classical phase space and the 
freedom in the (minimum) uncertainties in the canonical variables. We propose to consider 
states based on Voronoi complexes as playing a similar role for the discrete-to-continuum 
transition. We do not know yet how to phrase a minimum uncertainty condition in this 
context, or questions about the existence of processes which might produce such states. 
Filling in the first gap is one of the goals of our program; the second one will probably 
require a much more complete knowledge of semiclassical quantum gravity, including its 
dynamics. From a geometrical point of view, however, there is a strong motivation for using 
Voronoi complexes, that we will be exploring in this paper. 

Second, it is known that in order for states based on a set of graphs to span a dense 
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subset of the (kinematical) Hilbert space of loop quantum gravity, one needs to consider 
graphs with an arbitrary number of edges and connectivities. If one restricts oneself to 
graphs with only four-valent vertices (in 3+1 dimensions), one does not obtain but a high- 
codimension subspace of the Hilbert space. We are nevertheless suggesting that by restricting 
our attention to states defined over such graphs we will not lose important information. Are 
the semiclassical states defined over such a restricted class of graphs sufficient to display the 
needed semi-classical features? We do not have a definite proof for this, but we can argue 
in favor of such states. Consider for instance the example of a simple harmonic oscillator 
with a finite number of degrees of freedom. In this case the usual Gaussian coherent states 
one defines, peaked at phase space points, span only a finite-dimensional submanifold of the 
infinite-dimensional Hilbert Space of the theory. In spite of this, one can regard the coherent 
states as 'enough' for describing semiclassical states in some cases. Again for the coherent 
states of the free Maxwell theory, one can take coherent states and they approximate very 
well the semiclassical properties that we are interested in. We will then by analogy assume 
that the states we are considering here, defined over Voronoi complexes, will be enough to 
describe the semiclassical sector of the theory. 



III. THE INTRINSIC GEOMETRY OF VORONOI COMPLEXES 

Having introduced the necessary background concepts, we can now describe in more detail 
what we intend to do in the rest of the paper. Our general goal is the following: given 
a Voronoi cell complex Q, i.e., one satisfying the incidence relations ()2.2|) . determine the 
range of classical geometries that are consistent with it. From the quantum gravity point of 
view, it is not clear whether it is reasonable to associate a single fl with a semiclassical \1/ 
(different points of view underly for example the proposal in Ref 3|, or the shadow state 
proposal [15[, in which different discrete structures are seen as tools for probing the state 
We do so here because it allows us to separate the effects of the classical uncertainty 
in the discrete-to-continuum transition from those of the quantum fluctuations in and if 
more f2's need to be considered one can always combine their uncertainties later. 

We emphasize that the effective geometry will only be a spatial one; the recovery of an 
effective spacetime geometry requires either structures that can be interpreted as discrete 
spacetimes and the use of Lorentzian statistical geometry (see, e.g., Ref lii]), or additional 
variables on Q that can be interpreted as discrete versions of dynamical data (as in Ref [Sj). 
Even in this context, our actual calculations will concern cases in which the topology of Ai 
is trivially determined by that of Q, and only Sec II VI will discuss a more general situation. 



A. The Discrete-Continuum Transition 

There is an analogy with our situation in elementary physics. When one is dealing with a 
fluid, one 'knows' that at some 'microscopic scale', one is dealing with molecules, individual 
entities with which one can associate, classically, a position and a velocity. One then con- 
siders cells at a mesoscopic (crossover) scale inside of which one averages velocities, energies 
and so on, and one assigns such quantities to the cell as a whole. Finally, one goes to much 
larger, macroscopic scales and regards those properties of the cell as being local, deflned by 
continuum (and differentiable) flelds. In a sense, this is the procedure we are envisaging: 
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There is a microscopic (discrete scale) £d where the 'true' discrete geometry is defined by a 
graph. We will then consider large sets of cells at a scale £c, over which we will average the 
combinatorial quantities of the cells, to smooth out statistical fiuctuations and define meso- 
scopic quantities that vary slowly between such groups of cells. On the larger macroscopic 
scale, this will allow us to view the mesoscopic quantities as the local values of continuum 
fields that define a geometry. This passage between the microscopic scale C.^ and the final 
macroscopic one is what we call the discrete- continuum transition. 

To define a macroscopic geometry the quantities we consider will be geometric 

invariants Qx associated with extended submanifolds X (Z M., large enough to correspond 
to a large number of cells of VL considered as a cell decomposition of M. , but macroscopically 
small so that they do not correspond to integrating or averaging over regions where the geom- 
etry varies. In this paper, the submanifolds X will be simply open regons of }A containing 
large numbers of D-cells in D dimensions, although for other purposes jq] one might consider 
hypersurfaces in Ai, approximated by large collections of {D — l)-cells, or submanifolds of 
higher codimension. As for the Qx themselves, the above analogy with the thermodynamic 
limit leads us to divide the possible invariants into extensive ones (the simplest example 
is the volume Vx-, for which we obtain values by counting either Voronoi cells or Voronoi 
vertices contained in X) and intensive ones (examples of this type are curvature invariants, 
averaged over X). For the purposes of this paper, the latter are the less trivial and the more 
fundamental ones (relationships between most quantities are curvature-dependent), and we 
will concentrate on those in this paper. By constructing an effective continuum geometry 
here we thus mean finding the values of a sufficiently large set of curvature invariants Qx-, 
and providing a quantitative measure of the goodness of the construction, in terms of the 
values of the uncertainties AQx, using VL. 

In order to learn how to do this, we will first provisionally assume that a classical ge- 
ometry with known values for all of its Qx is given, which allows us to calculate statistical 
distributions of Voronoi complex variables. The relationship will then be inverted to allow 
us to estimate the Qx and their uncertainties from a given Vt. In other words the main idea 
is that, when one obtains a cell complex as the result of applying the Voronoi construction to 
a random point sprinkling of density p in a Riemannian manifold (A^,5'afe), the combinato- 
rial quantities for the complex satisfy dimension- dependent identities, and one can also 
calculate (at least in principle) geometry-dependent probability distributions for values of 
those quantities. These two types of relationships together imply that the complex encodes 
enough information about the manifold M. and metric qah that we could reconstruct (A^, qab) 
from it, up to statistical uncertainties on volume scales at or below p~^. 

Thus, if we are given an abstract, embeddable Voronoi complex in the sense of Sec m 
we can construct an approximate (Al,^^;,) that is a good continuum version of Vt on scales 
larger than the average embedded cell size VmI^d ='■ P~^- From a practical point of view, 
we quickly run into the difficulty that only a few of the relevant probability distributions, 
for low- dimensional fiat or constant-curvature spaces, are known. Therefore, we will treat in 
detail the two-dimensional case, where we can derive the results we need analytically, and 
outline the procedure in three dimensions, where analogous calculations will have to be done 
with computer simulations. 
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B. The Two-Dimensional Case 



In two dimensions, any metric is conformally flat, and can be locally written as qab = Gab, 
where / is a scalar function, and Cab a fixed flat metric on a portion of the 2-manifold; if 
Cartesian coordinates are used for Cab, the line element is then qab da;"da;'' = e^^^ (dx^ + dy^). 
The function / is in turn related to the scalar curvature hj R = — 2e~'^V^/; R is the 
continuum geometrical quantity we will associate with a cell complex in this section. When 
attempting to construct a geometry, it would seem natural to consider first finding a distance 
for any pair of cells (or of vertices) in the complex. However, the relationship between 
areas and lengths is curvature-dependent, and if we consider the cell density to be a basic 
parameter, before we assign distances to pairs of objects in the cell complex, we first need to 
find out the curvature that best fits each portion of Q. Besides, in physical applications one is 
often directly interested in the curvature of a manifold, since it affects, e.g., the propagation 
of matter fields on it. 

In order to find out how to assign a value of i? to a subset of a cell complex Q, we need 
to learn to recognize complexes that might arise from a point sprinkling in a manifold with 
curvature R. We will therefore use Eqs (12.11) and ()2.3|) to determine the mean and variance 
of the number of edges of a 2-cell on a 2-sphere (S^, Sab) of constant scalar curvature. If we 
denote by Qr a Voronoi complex on such an S^, Eq ()2.1|) becomes A^o — + ^2 = x(^-r); 
then, using ()2.3|) and noticing that the average value of A'^i(ci;) over the complex is given by 
{Ni{u))^^nj^ = 2 Ni{Qji) / N2{^r) , where the 2 is due to the fact that each edge is shared by 
two 2- cells, we obtain 

^)-(- a- 

Here, we have used the Gauss-Bonnet theorem relating x{^r) and the scalar curvature of 
(S^, Sab)- Since the average does not depend on Qr, it also equals the mean number Ni of 
neighbors of a cell taken over all random complexes in the geometry (S^,Sab) with density 
p. Eq (j3.1|) can then be easily inverted to give the scalar curvature in terms of the mean 
number of edges of a cell, 

i? = 47rp(l-|iVr) . (3.2) 

Finding the variance af = Ni — (A^i)^ of the distribution of the number of neighbors is 
considerably more difficult. We already know (A^i)^; for Ni, we will use a trick Notice 
that for any 2-cell uj, Ni{uj) equals the number of vertices Nq{u), so our task can be seen 
as that of calculating N^. The latter has been calculated (i.e., analytically reduced to an 
integral which is then numerically evaluated) by Brakke [l7[ for the fiat case; we will now 
generalize his calculation to the case of a constant positive curvature manifold. For any uj, 
a quantity related to whose mean is much simpler to calculate directly is the number of 
(unordered) pairs of vertices not sharing an edge, for if we denote this quantity by A^o,o'('^)5 
then by definition A^'q^o' = {^o ~ 3). Thus, Nq = 2 A^o,o' + 3 iVo (notice that our Nofii is 
Brakke's I{v,v)), where Nq^/ can be found by integrating a suitable probability density, as 
we now show. 

Consider the 2-sphere as embedded in 3-dimensional Euclidean space, and call C the 
center of the sphere; the radius a of the sphere is related to the scalar curvature of by 
R = 2/a^. Choose an arbitrary point O on the 2-sphere as its origin. To locate any other 
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point P in we will initially use the spherical coordinates (Xi^)) where x ^ [0)'^] is the 
angle at C between the lines CO and CP, and 9 G [0, 27r] the azimuthal angle on around 
O. The line element on is then given by the familiar form 

ds^ = Sab dxMx'' = (dx^ + sin^ x d^^) . (3.3) 

Consider now an arbitrary cell in a random Voronoi complex on (S^, Sah) of density p, and for 
convenience choose the coordinates such that the sprinkled point or "seed" 5*0 that defines 
this cell is at the origin. We would like to find the expected number of pairs of vertices 
(Pi,P2) of this cell which do not share an edge, i.e., which are not consecutive. Any cell 
vertex Pj is equidistant from three seeds in the sprinkling, in this case Sq and two others, 
{Sii,Si2), where i ~ 1,2 labels the vertex they are associated with. Therefore, when we 
count pairs of vertices we need to count pairs of configurations (5*0, 5*^1, 6*42) such that the 
disk Vi inside the circle through each triple of points is void of other seeds (so that they 
really define a vertex) , and the four seeds Sij are distinct (so that the vertices Pi and P2 are 
not consecutive). What is the probabihty density for all of this to happen, in terms of all 
possible locations for the four seeds in question? 

The probability measure for a seed to be located at S* = (x, 9) is p (a^ sin x dx d^) = 
pa^ d(cosx) d6'; therefore, for the pair of seeds {Sii,Si2) giving the vertex Pj, at locations 
iXii, 9ii) and {xi2, 6*12), respectively, with i = 1 or 2, 

= p^a^ d(cos xn) d6'a d(cos Xi2) <^9i2 . (3.4) 

Wc can make sure that none of the four seeds is contained in the disk defined by and the 
two seeds in the other pair by specifying appropriate ranges for the allowed positions of the 
seeds, and we can impose that the union of the two disks contain no additional seeds by 
multiplying the measure by the probability Q-P-^i^'^^^'^) , 

To control more easily the ranges of integration, it is convenient to make a coordinate 
transformation from the eight variables (xij, %) to a set of four variables (C, 0, t^i, 002) which 
specify the location of Pi and P2 (or equivalently, the size and location of the two circles), 
and four variables aij which specify the location of the four points on the two circles. Given 
two circles through Sq, or two points P, = (x?, 9i) , labeled so that 62 — 9i < tt, call Q — (2^, 0) 
the other point at which the circles intersect, besides 5*0 (if the circles are tangent at we 
identify (j) with the direction of that tangent, and Q = So, but this happens with probability 
zero). Also, call Ui and UJ2 the angles QSqPi and QS0P2, taken to be positive respectively 
in the clockwise and counterclockwise directions from 0. Then 

9i = (I) - uji , 92 = (I) + UJ2 , (3.5) 
and the ranges of values for the new angles are 

< (/)< 27r , -| < a;i < ^ , -uoi < LJ2 < ^ . (3.6) 

The distance Xi of Pi from Sq can be expressed in terms of ( using the spherical cosine law 
applied to the isosceles triangle SqQP^, i.e., cosxi = cosxi cos 2^ + sinxi sin2C cosa;^, which 
implies the relationship 

tanxi = ^, with < C < 77 • (3.7) 

cos uJi 2 
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FIG. 2: Construction for the calculation of 2D statistical fluctuations on the two sphere; see text 
for explanation (only one of the two Pi and one of the Sij are labelled) . 

The additional variables are defined as follows. On the circle around Pi (resp, P2) call 
5*^1 the first seed one meets going around the circle clockwise (resp, counterclockwise) from 
5*0, and Si2 the second one. Then an is the angle PiSoSn and the angle PjS'o5'j2 (both 
measured clockwise for i = 1 and counterclockwise for i = 2), i.e., 

% = 9, + i-iyaij = + (-1)^ {coi + aij) , (3.8) 

with ranges given by 

TT TT 

~uji< an < - , an < ai2 < - . (3.9) 

The Eqs ()3.8p constitute half of the coordinate transformation. To find the other half, we 
apply the spherical cosine law to the triangle SoPiSij] this gives cosxi = cos Xi cos Xij + 
sin Xi sin Xij cos aij, which implies the relationship 

tan ^ = tanxi cos aij . (3.10) 

We can then substitute ()3.7|1 in this last equation, and express the result in the more conve- 
nient form 

cos^ uji — tan^ ( cos^ a^,- zo -, 1 \ 

cosXij = — ^ , , 2 • (3-11) 

cos^ uJi + tan (, cos^ aij 
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With ()3.8|) and ()3.1H) we can now rewrite the full measure of integration ()3.4p in terms of 
the new variables. After a somewhat lengthy calculation, one obtains 



d/ii d/X2 = dC d0 dui duj2 dan dau da2i da22 
X 256 tan^C(l + tan2C) 



n 



COS a,- 



cos^ Ui + tan^ C cos^ ^jj)^ 



X cos^cui cos^a;2 sin(ci;i + 102) sin(a;22 — tt2i) sin(a;i2 — an) • 



(3.12) 



We now need to calculate the area AiVi U V2) of the union of the disks around Pi and P2. 
The line SqQ divides each of the Vi into two parts, and Vi U V2 is the disjoint union of one 
part from each Vi (see figure) which, for the purpose of finding its area, it is convenient to 
think of as Vi with the "wedge" S^PiQ removed and replaced by the triangle S^PiQ. Thus 
we can write 

AiVi U V2) = Xli=i(^disk - ^wedge + Ariangle)i = Hi=i + 2u; - (vT - /?) COS x] , (3.13) 

where some simple spherical geometry gives 



^disk i 
-^wedge i 



A 



triangle i 



27r (1 — COS Xi) 0^ 

= 2(3i (1 - cosxi) 

2{(3i + uj^-\ vr) . 



(3.14) 



Here, (3i is half of the internal angle of the wedge or the triangle S^PiQ at Pj. Using the sine 
law with half the triangle QSoPi, we obtain that it is related to our variables by 



cos Pi = cos( sinwj , 



and we obtain from ()3.7j) that 



cos Xi 



COS^ UJi 



cos^ Ui + tan^ ( 
Putting these pieces together we therefore have 



(3.15) 



(3.16) 



A{Vi U V2) 



E 



n + 2u!i — 



2 cos uji 



COS^ LUi 



tan^ ( 



(vr — arccos(cosC sincijj)) 



(3.17) 



and the expectation value we are looking for is obtained integrating the product of these 
probabilities over all locations of the four seeds. 



■pA{ViUV2) 



Finally, substituting this in the expression for the variance of A'^o, we get 



af = 2 A^o,o' + 3 iVi - iVi = 2 / d/ii / dfi2 e 



-pA{ViUV2) 



- 18 1 



3 R ^ 1 
4:71 p Svr^ p2 



(3.18) 



(3.19) 
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This integral can be evaluated numerically for any given value of p/R. 

We are now in a position to discuss the geometries we associate with a 2-dimensional 
Voronoi complex f2. Since in two dimensions the curvature is completely characterized by 
the Ricci scalar i?, our goal is to associate a value of i?, with suitable uncertainties, with 
every set ?7 in a cover of a manifold 7W ~ i7 that is made of sufficiently small sets. We 
start by selecting a larger collection of candidate sets, and then explain how to pick the 
appropriate ones among them. For each 2-cell ujq E VL consider the family of sets {U^^^^^\}, 
A = 0, 1, 2, defined by 

= ^0, U^,,x+i = U^o,x lJ{ioEn\dujn dU^,,x 7^ 0} ; (3.20) 

in other words, U^q^x is uq together with the first A layers of neighboring 2-cells around it. 
Roughly speaking, if we call A^a;o,A the number of 2-cells in U^^o^x, in an approximately fiat 
2-geometry we expect to have N^^q^x ~ 1 + 6 + 12 + ... + 6A = 3A^ + 3A + 1. A value of R can 
be assigned to each c^o simply by using ()3.2p as an estimate, 

R{uJo)=Anp{l-lN,{LUo)). (3.21) 

This value however is not a good one as far as manifold geometry goes; differences between 
it and those obtained for neighboring 2-cells should be interpreted not as real variations of R 
but as statistical fluctuations. What we should do instead is average the values obtained for 
a cluster of neighboring 2-cells, i.e., a suitably large Ui^o,x- If we knew that the continuum 
geometry will turn out to be a constant curvature one, the best strategy would be to pick 
A as large as possible, in order to maximize the statistics. In practice, if we use too many 
2-cells we may be combining regions that in a good flt of the geometry would have different 
curvatures, so we need to specify a procedure for picking an optimal A. Pick some value of 
A. Then, at scale A, ?7^o,a provides a sample of size N^^^x from the ensemble of all 2-cells 
in a manifold of assumed constant curvature, and our best estimate for the scalar curvature 
around ujq is the average 

R{U^o,x) = = 47rp (l - 1 {N,{uj))^^u^^^^) , (3.22) 

with the variance in the distribution of averages over such samples being given by 

^W) = ^Jo!a^?- (3.23) 

If this variance is small, we can use it to estimate a range of values of R which would have 
produced values of (A^i(co')) within this tolerance, by setting AR = {dR/d{Ni)) cr(Arj), or 

Ai?(t/.o,A) = ^-(^.) = ^<f-i- (3.24) 

Finally, the whole process is self-consistent and gives a good continuum approximation to 
Q if differences between the scalar curvatures estimated from neighboring regions U^q^x are 
smaller than the statistical uncertainty AR{Uujg,x) within individual regions, i.e., R is slowly 
varying on the scales given by A. 
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C. Considerations on the General Case 

In three dimensions, there are three independent local curvature invariants. Different choices 
are possible; for example, they can be thought of as the eigenvalues of the Ricci tensor Rat, 
or as the scalars R = R^, R%R\, and R%R\R^a- To recover them from VL we need to use 
three independent local combinatorial graph quantities, chosen among the Nk{u). One may 
think that a possible choice for these is the set of numbers {No{C) , Ni{C) , N2{C)) , which 
count the number of vertices, edges and 2-faces of a 3-cell C. However, for any given C, 
only one of those three numbers is independent, since its vertices, edges and faces form a 
2-dimensional complex dC homeomorphic to S^, to which we can apply Eqs I2.1H2.3I with 
x{dC) = 2, and we get that iVo(C) - Ni{C) + N2{C) = 2 and Ni{C) = lNo{C), which is 
easily solved to give, for example, 

No{C) = 2 N2{C) - 4 , Ni{C) = 3 N2{C) - 6 . (3.25) 

It seems reasonable to assume that, analogously to the 2D result, the number N2{C) of 
2-faces of a 3-cell allows us to estimate the scalar curvature R. The other curvature invari- 
ants are associated with the anisotropics in the curvature, since knowledge of the traces of 
powers of i?"^ can be thought of as corresponding to knowledge of its eigenvalues. We can 
then conjecture that this type of curvature information is encoded in direction-dependent 
combinatorial quantities of the 3-cells, i.e., any statistically significant differences between 
properties of 2-cells on different parts of the 3-cell boundaries; for example, the direction 
in which faces of a 3-cell C have the least number of edges are those in which dC is "most 
curved", and should correspond to an eigenvector of R'^b- Checking such statements is be- 
yond the scope of this article, and will most likely require computational work with families 
of simple 3-geometries. 

IV. SMALL-SCALE GEOMETRY AND COARSE-GRAINING 

In the previous section we assumed that all complexes are embeddable in some manifold, 
but one of the ultimate goals of graph-based approaches to quantum gravity is to obtain 
a manifold-independent formulation. This section represents a tentative first step in that 
direction. It is somewhat more speculative than the previous ones, both for this general 
reason and for a more specific one. Even though our motivation comes from the semiclassical 
sector of quantum gravity, we will continue to work in a purely geometrical setting, i.e., we 
do not consider any quantum aspects, states, operators and the like. What we shall assume is 
that it is meaningful to consider a fixed discrete structures even at scales below the cross-over 
scale we have considered before, either because the quantum-to-classical transition occurs at 
a still smaller scale £q < l^i oi because we can anyway learn something about the continuum 
interpretation of the theory by pushing this geometrical view as far as possible. This is the 
meaning we shall give to the term "small scale geometry" . We will discuss two aspects of this 
issue here, the different scales ("combinatorial" scales C in the discrete structure, or metric 
scales L in the continuum) that play a role in the construction of a Riemannian manifold, 
and how to deal with some of the obstructions that may come up in this process. As we will 
see, even in the simple setting in which we consider only the milder obstructions, different 
situations may arise; we will need to introduce the concept of coarse-graining, and to define 
other length scales in addition to the ones we have already seen in other sections. 
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A. Length Scales 

One general issue that needs to be addressed is that of identifying the various length scales 
that appear in our framework, and assigning dimensional values to combinatorial quantities 
in Q. Since we are assuming that the fundamental theory is formulated in terms of Q, we 
would want to describe those length scales using quantities intrinsic to Q, and be able to 
make statements such as "at scale £q, Q cannot be embedded in a manifold," which implies 
that a continuum geometry is not available to provide a meaning for length scales. 

Continuum theories include a number of length scales associated with various phenomena; 
results of measurements are interpreted in terms of those scales and provide values for them, 
so in order to provide values for scales intrinsic to Q it is natural to go through a correspon- 
dence with continuum length scales. An obvious one to try to identify in terms of Q is the 
Planck length Lp (when Q is dressed with other variables it may be possible to do this even 
at the kinematical level, as in the case of the calculation of the loop spacing for the heuristic 
weave states in loop quantum gravity but any prediction of the discrete theory that 

has a continuum counterpart can be used, and when more than one will be available the 
discrete theory can be tested. 

A discrete structure of the type we are discussing also provides various scales. For example, 
if we start with a Voronoi cell complex, we can embed it in a geometry that has no length 
scales smaller than the cell size. The latter cannot be assigned a dimensional value yet, 
because the situation is invariant under a global rescaling of the metric, but we will say 
that Q is associated with a semiclassical cross-over scale 4 equal to its discreteness scale 
id- If we start with a D-dimensional Voronoi complex Q that is not a cell complex, but 
is such that we can obtain a cell complex Q' by coarse-graining Q and removing a fraction 
< ^ < 1 of its D-dimensional tiles, we will say that Q' corresponds to the semiclassical 
scale £s, at which it and its dual Q'* represent the manifold well, while Q is characterized by 
the discreteness scale £d = (1 — O^^^^s- In either case, it may not be possible to determine 
the metric or curvature of the geometry at scale is, because at that scale we do not have 
enough statistics to reliably use the techniques we will describe. Instead, the classical scale 
ic is the linear size, in terms of number of cells, of the regions in Q used as approximations 
to the manifold regions X such that the continuum geometry can be determined with small 
statistical fluctuations. That geometry itself, possibly together with the global topology of 
M., may also determine larger length scales ig (both local ones defined by the curvature 
and its rate of change or higher derivatives, and possibly global ones defined for example by 
non-trivial homotopy generators). 

Summarizing, the length scales < 4 < 4 < ig are determined by the structure of Q, if 
$7 is a discrete structure of a type that allows us to define 4, given a suitable definition of a 
coarse-graining procedure. In quantum gravity, we will then require, as part of the definition 
of a semiclassical state, the condition that be associated with a discrete structure for 
which 4 exists. On the other hand, continuum-based physics determines a length scale 
Lp for quantum gravity, given in terms of fundamental physical constants by dimensional 
arguments and back-of-the-envelope calculations of quantum gravity fluctuations 0, which 
appears also in the eigenvalues of the quantum geometry operators for areas and volumes; It 
is not clear to us whether this scale should be identified with 4 or 4, if different, or how this 
relationship might be affected by renormalization arguments. Finally, physics also determines 
a phenomenological, "macroscopic" scale ^ Lp, which depends on the experimental 
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techniques used but that can be assumed to be larger than 4- A good understanding of the 
relationships between all of these scales is outside the scope of this paper, but our discussion 
of coarse-graining will be a start in this direction. 

B. Coarse-Graining 

Some complexes Q can be directly associated with a "macroscopic" geometry, in the sense 
that they can be embedded homeomorphically in a manifold at a length scale 4 and the 
procedures previously described give a well-defined (approximate) geometry at a length scale 
is, in which case the "semiclassical" interpretation is relatively straightforward, within the 
appropriate statistical uncertainties. That is, Q gives us a classical, mesoscopic geometry. 
At larger scales, for which the number of points to be taken per sample cell is larger, one 
expects a very good approximation, according to our own set of conditions. The remaining 
question then is: What if fl does not satisfy the requirements for approximating a nice 
geometry? Although there then is no direct association between graph and geometry, if 
in an appropriate sense the obstructions only occur on small scales, it may be possible to 
coarse-grain the graph to a larger-scale one that is associated with a continuum geometry. 

The coarse-graining procedure will be one which, intuitively, takes a graph at a certain 
scale and maps it to a different, larger-scale Q'. There are two possibilities on how to perform 
the coarse graining in this case. The first possibility, that we call soft coarse-graining, has 
the following strategy: One assumes that one has a Voronoi or Delaunay complex and one 
gets to a coarser one by "ignoring" the structure below a certain scale simply by averaging 
the quantities of interest over neighboring cells. This is what we have already done in Sec lIIII 
when looking for the optimal mesoscopic scale at which to define slowly varying quantities; 
we are not actually defining a new complex, and it can always be done. In the second 
possibility, that we will refer to as hard coarse-graining, we have a procedure that connects 
two graphs, and one gets the larger graph Q' by removing cells (Delaunay vertices) from Q, 
in a precise way. 

Hard coarse graining, which implies a modification of the underlying graph, represents a 
new input in our considerations. This procedure is achieved by a sequence of cellular moves 
that refine the cellular decomposition [l9|. These moves have been shown to provide a pre- 
scription for refining and coarsening the complex (they are in a sense more elementary than 
the Pachner moves). The intuitive idea for the coarse graining is that two adjacent n-cells of 
the complex that share an [n — l)-cell get 'fused' into a new n-cell. This refining/coarsening 
procedure has already been used for defining renormalization prescriptions in discrete sys- 
tems P^, and the approximation to smooth manifolds by our proposed hard coarse-graining 
is an explicit example of some of those prescriptions. 

In order to have a proper understanding of the different scenarios that one might en- 
counter given an arbitrary graph at small scales, with the purpose of incorporating it into 
our geometrical formalism, let us make the following classification in decreasing order of 
'complexity': 

1. In the most exotic scenario one starts with a Voronoi complex that does not have a dual 
cell complex. This means that there is no dual triangulation and the corresponding 
graph might be disconnected. Since this happens for cases where, say, one sprinkles 
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points in a manifold with a density that is "too low" , and no real statistics is possible, 
we will exclude this case from our considerations here. 

2. The next possibility is that the given cell complex Q can not be embedded into any 
manifold A4. This means that the complex is not topologically a cell decomposition of 
a manifold. In this case one might hope that with a proper coarse graining one might 
be able to take this complex to a new one Q' that does admit such interpretation and 
can be embedded. Physically, one could say that the original complex represented a 
space that on certain small scales has 'topological' and/or 'dimensional' fluctuations. 
This is the trickiest situation where one needs a 'hard coarse graining', for which we 
are not aware of an existing specific procedure, and it is certainly a question worth 
pursuing. 

3. A milder scenario is to consider as the starting point a cell complex Q that is a cell 
decomposition of a manifold A4 but that is not Voronoi, in the sense that the valence 
of the graph is not equal to D + 1. In this scenario, one would allow for "fluctuating 
valence" in the Voronoi cell, but with average (A'^iio) = 4, in 3D. This means that 
the graph that defines the quantum state at Planck scale does not have the 'Voronoi 
signature', which might be the case for an arbitrary graph in loop quantum gravity. 
In this case we can implement the cell moves of Ref |19] to turn it into a Voronoi 
cell complex; if the authors' conjecture is correct, via these cellular moves one can 
always transform any cell decomposition into a Voronoi one, that will serve as starting 
point for the statistical geometry considerations described in previous sections. One 
should also note that in some cases this 'restructuring' of the cell complex might not 
correspond to a course graining in the strict sense, since the moves that 'lower' the 
valence of the Voronoi graph do not remove cell complexes to it as one might expect 
from a coarse graining procedure. 

In the next, final section we shall consider quantum issues and relate the formalisms we 
have been employing to the task of bridging the gap between a semiclassical state in loop 
quantum gravity and classical geometry. 



V. PLANCK SCALE GEOMETRY AND BEYOND 

In this section we make precise the connection between the previous sections and semiclas- 
sical canonical quantum gravity, i.e., the low energy limit of loop quantum gravity. This 
section has two parts. In the first one, we outline how one can try to use our framework to 
complement present approaches to semi-classical states in LQG. In the second part we outline 
the steps to be followed in order to estimate statistical errors in making phenomenological 
predictions out of the semiclassical states. 

A. Semiclassical States 

In order to deal with the low energy limit of the theory, we will again have to deal with the 
fact that there might be new length scales in the description of the quantum geometry. So 
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far we have considered in great detail the way in which we may approximate a macroscopic 
geometry starting from a graph or complex at a much smaller scale. Even though we have 
introduced a microscopic discrete scale I a-, at which the finest graph might be defined, this 
scale is not necessarily related to a particular quantum scale £q. As discussed in SecHH there 
are standard dimensional arguments which indicate that the relevant scale £q for a quantum 
state is the Planck length Lp. This only means that the quantum excitations of the geometry, 
which could be elements of area through edges of the graph or contributions to the volume 
from the vertices, yield eigenvalues of the order of the Planck scale for single edges/vertices 
(at least for the unrenormalized value of the Barbero-Immirzi parameter 7 of order one that 
is used nowadays). Whenever we consider a semiclassical state \l/sc in loop quantum gravity, 
there are several steps and assumptions involved in estimating the continuum geometry it 
corresponds to. Let us now explore what those assumptions are. 

The first assumption is that the true nature of geometry at the Planck length is described 
by loop quantum geometry. Second, we assume that we can define certain states (Coherent 
2^ . Shadow Gaussian ||26j, etc.) that will be semiclassical in the sense of approximating 



a geometry (and its time derivative), in the 'correct phase'. Finally, one is assuming that 
a definite picture of the quantum geometry will emerge from a precise merger of both loop 
quantum gravity and our "discrete-statistical" approach. In this viewpoint, the final picture 
looks like this: The 'true state' is given by a properly defined state in loop quantum gravity, 
featuring quantum behavior at the Planck scale. The graph that we use to define this 'shadow 
state' (to give an example) will be a Voronoi graph. The semiclassical state will have the 
property that, when probed on a mesoscopic scale 4 (still to be specified in practice) it 
will behave quasi-classically (hopefully, with respect the the coarse-grained operators that 
we need to define). This means that in order to construct the desired semiclassical state 
following the steps of Ref [2l[, the operators that we specify as belonging to the set to be 
approximated (together with the tolerances) will not be operators assigned to observables 
at the quantum, Planck scale. 

We are interested in constructing states that approximate quantities at a scale 4, where 
the quantum-classical transition takes place. As already mentioned, this new scale, yet to be 
identified, could be assumed to be close to the mesoscopic scale 4, but it should be smaller. 
Going to our analogy model of fluids, the transition quantum-classical can be assumed to 
take place at a smaller scale than the one at which averages are taken and the continuum 
approximation emerges (for the simple reason that one adds velocities of particles contained 
in the region, which already presupposes classical attributes for the constituents of the fluid). 
The viewpoint that new, mesoscopic scales have to be considered has already been explored 
by other authors who have explored the semiclassical limit before [2^. 

The next question we would like to address is the possible utility of our statistical methods 
for distinguishing between good and bad candidates for semiclassical states. That is, suppose 
that we are given an alleged semiclassical state \l/sc and we are assigned the task of testing 
it. The strategy for doing so is the following. Compute expectation values and fluctuations 
for the cell-based, observables that the maker of the state claims it approximates well. Make 
sure that the state is defined on a Voronoi complex. Dress the complex with the expectation 
values found for the observables. Now the crucial step is to compare this dressing with the 
geometry we expect to get from the complex itself, as in Sec IIIH and verify that they are 
consistent. If they are not, return the state to the maker. 
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B. Phenomenology 



Suppose now that the given state satisfies all criteria set by some considerations, then in 
our procedure we would have to pass it to the next step in the testing line: see whether it 
fits observations. Now this is a very intricate question, for we have not specified the phe- 
nomenological criteria that a semiclassical state should satisfy, nor the expected observable 
deviations from the classical realm (Lorentz violation, modified dispersion relations, non- 
commutativity, birefringence, etc) 122]. Note also that this procedure implies the knowledge 
of more quantities assigned to the cell complex, more than what we have up to now specified 
in the complex. This further 'dressing of the complex' will be dealt with in a forthcoming 
paper 13]. 

What we could attempt to answer is the following simpler question. Suppose that we 
expect the state to approximate the geometry gat- As we have discussed before, there are 
now two sources of uncertainty in the approximation of this spacetime geometry (even though 
we were working in the canonical picture, we can talk about a spacetime geometry). The 
first source is a well known one, due to quantum fiuctuations: observables 0/ that were 
chosen to be approximated by the state have quantum dispersions (AO/)^. The second 
source of uncertainty is the statistical nature of our reconstruction procedure. That is, as 
we have argued before, the macroscopic continuous geometry is only a fiction, very useful 
for describing phenomena at certain scales, but it is only an approximation to the true 
geometry. But there is not a unique macroscopic metric that can be approximated by the 
graph: rather, there is an 'ensemble' of such metrics, which then introduces a classical 
or cross-over uncertainty (AOj-)c in the observables. The observational imprint that this 
uncertainty might bring is manifested by the fact that we normally use a fiducial classical 
metric to perform geometrical measurements. We expect that this geometry belongs to the 
set of possible geometries that are approximated by the graph, but it might turn out to 
be not the 'most probable' one, and thus one might be introducing an extra error in our 
interpretation of the experimental results. Note that these errors are of a different nature 
than those that arise from the use of a regular lattice in We expect that we won't have 
these 'discretization errors' due to the use of the Voronoi construction. In order to have 
full control over the different sources of uncertainty, one would like to quantify them and 
decide which one is the dominant one. At this point we are not in the position of having a 
working hypothesis for this question, and we can only make an educated guess as to the way 
in which the two uncertainties interact with each other. Under the most naive assumption 
of independence, one might say that the total uncertainty is 

{AO,)l = (A0,)| + {AOj)l . 

An important issue is to quantify the fluctuations that can be attributed to both effects, and 
know how the total fluctuation for each observable depends on them, since we would like to 
be able to subtract any statistical 'spurious' information and be able to measure the pure 
quantum contributions to the problem. 

Assuming that the semiclassical state had these desired properties — an issue that lies 
outside the scope of this paper — would lead us to conclude, from a strict viewpoint that 
we have succeeded in our task of constructing a semiclassical geometry. This is because 
this state should then approximate any observations, at the mesoscopic level that we shall 
perform on the geometry (with external matter) and thus, if we now take the viewpoint that 
it is enough to describe what we can observe, we would reach our conclusion. 
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C. Summary and Outlook 



To summarize, we have introduced a framework for establishing a bridge between the de- 
scription of geometry of space in terms of semiclassical states for quantum gravity, and the 
one used in phenomenological calculations of corrections to field theory in curved space. At 
the most basic level, our proposal fits in with the general idea that the type of discreteness 
one encounters in a class of approaches to quantum gravity implies that the amount of infor- 
mation contained in a finite spacetime volume is finite, which has been explored in a variety 
of contexts . Our work then suggests a specific way to translate this idea into quantitative 



expressions for the information a discrete structure contains about a continuum geometry, 
using statistical geometry. To illustrate how the framework is used, we have characterized 
a semiclassical quantum gravity state by a single cell complex and a single set of values 
for variables motivated by loop quantum gravity. Such a characterization, of course, is not 
intended to be complete; it might arise as part of the specification of a coherent state f24l |. 
or from a "shadow state" used to probe "ifgc but for our purposes it is just a tool that 
allows us to isolate the statistical effects we want to study. It is a tool that complements the 
criteria one might have in choosing the quantum state. 

The general picture that emerged in this framework is one in which, from a certain length 
scale upward, if the quantities characterizing the state are sufficiently regular, the mesoscopic 
geometry is determined up to quantifiable uncertainties. Once those uncertainties will be 
better understood, especially in three dimensions, one can begin to look into the possibility 
of adding matter coupled to gravity, and how these uncertainties need to be combined with 
quantum fluctuations to give measurable effects on observable quantities. The answer to this 
last question will probably depend on which observations are carried out, and it is in principle 
possible that one can use observable effects that depend only on quantum uncertainties; 
after all, going back to the gas analogy, when we want to study the quantum properties of 
atoms, we do not need to combine their quantum fluctuations with the statistical mechanics 
description of a gas. But a gas is different, and we do not need to have information about 
the state of the whole gas to even be able to describe what each atom does. 

In this paper we have concentrated our attention to graphs that did not have any other 
geometrical information other than what is already in its definition. In order to have an ap- 
proach closer to the variables used in LQG one would like to have assigned certain (classical) 
quantities to the graph Q, that is, one would like to consider 'dressed' graphs. This will be 
dealt with in a subsequent publication (8|]. 
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